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Abstract 

r  '  v 

Wfr  study  single  server  queueing  models  where  in  addition  to  regular  arriving 
customers,  there  are  negative  arrivals.  A  negative  arrival  has  the  effect  of  removing  a 
customer  from  the  queue.  The  way  in  which  this  removal  is  specified  gives  rise  to 
several  different  models.  Unlike  the  standard  FIFO  GI/GI/1  model,  the  stability 
conditions  for  these  new  models  may  depend  upon  more  than  just  the  arrival  and 
service  rates;  the  entire  distributions  of  interarrival  and  service  times  may  be  involved.  ' 

rx 


Key  Words:  queue,  negative  arrival,  stability. 


August  20,  1989 


1Ecole  des  Hautes  Etudes  en  Informatique,  Universite  Rene  Descartes  (Paris  V),  45  rue 
des  Saints-Peres,  45006  Paris,  France.  Research  supported  in  part  by  CNRS-CS 
(French  National  Program  on  Parallel  Computation). 

department  of  Operations  Research,  Stanford  University,  Stanford,  CA  94305-4022. 
Research  supported  by  the  U.S.  Army  Research  Office  under  Contract  DAAL-03-88-K- 
0063. 

department  of  Industrial  Engineering  and  Operations  Research,  Columbia 
University,  NY,  NY  10027.  Research  supported  by  NSF  Grant  DDM  895-7825. 


1 


1.  Introduction 


Consider  a  single  server  FIFO  queue  with  two  types  of  arrivals,  regular  and  negative.  Regular  arrivals 
correspond  to  customers  who  upon  arrival,  join  the  queue  with  the  intention  of  getting  served  and  then 
leaving  the  system.  At  a  negative  arrival  epoch,  the  system  is  affected  if  and  only  if  customers  are  present  in 
which  case  a  customer  is  removed  from  the  system.  Intuitively,  the  introduction  of  negative  arrivals  makes 
the  system  less  congested  than  if  they  were  not  present.  In  particular,  one  would  expect  that  a  steady-state 
distribution  for  queue  length  can  exist  even  when  the  regular  arrival  rate  is  greater  than  the  service  rate. 
In  the  present  paper,  we  consider  several  single  server  models  of  this  type  and  give  necessary  and  sufficient 
conditions  for  the  existence  of  a  unique  steady-state  distribution  of  queue  length.  As  we  show,  the  stability 
conditions  may  depend  upon  more  than  just  arrival  and  service  rates. 

The  notion  of  negative  arrivals  was  introduced  in  [3]  in  the  case  of  a  Markovian  network  generalising 
a  Jackson  network  [4]  to  include  negative  as  well  as  regular  (positive)  arrivals.  In  particular,  it  was  shown 
[3]  that  such  a  network  has  product  form,  though  its  customer  flow  equations  are  non-linear.  Several 
practical  applications  motivate  these  kinds  of  models.  Negative  arrivals  can  represent  commands  to  delete 
some  transaction,  as  in  distributed  computer  systems  or  databases  [1]  in  which  certain  operations  become 
impossible  because  of  locking  of  data  or  because  of  inconsistency.  Negative  and  positive  customers  may  also 
represent  inhibitory  and  excitatory  signals,  respectively,  in  mathematical  models  of  neural  networks  [5,6], 
while  queue  length  in  this  case  represents  the  input  potential  to  a  neuron. 

2.  Removal  of  the  customer  in  service 

In  this  section  we  consider  the  case  when  a  negative  arrival  removes  the  customer  (if  any)  in  service  (RCS). 
We  will  study  two  different  models  of  this  kind.  The  first  model  assumes  that  regular  arrivals  form  a  marked 
point  process  tp  =  {(t„,Sn)  :  n  >  0},  where  ( t„,S„ )  denotes  the  n,k  customer’s  arrival  and  service  time. 
Negative  arrivals  form  a  Poisson  process  (at  rate  a)  assumed  independent  of  ip.  Let  L(t)  denote  the  number 
of  customers  in  the  system  at  time  t.  Let  {Znj  and  Z  be  i.i.d  ~  ezp(a)  independent  of  ip.  Our  first  model 
is  very  elementary  as  the  following  lemma  shows. 
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Lemma  2.1.  (L(t)}  has  the  same  distribution  as  if  the  model  was  a  FIFO  single  server  queue  with  marked 

point  process  ip  =  {( tn,Sn )  :  n  >  0}  with  5„  =  min(Sn,Z„). 

Proof:  Because  of  the  properties  of  the  Poisson  process,  when  the  nth  customer  begins  service,  the  amount 

of  time  spent  with  the  server  is  distributed  as  min(S„,Zn).  It  follows  that  {£({)}  is  stochastically  identical 
to  that  for  a  FIFO  single  server  queue  having  customer  arrival  times  {<„}  and  service  times  {S„}.  ■ 

The  importance  of  the  above  proposition  is  that  it  allows  us  to  analyze  our  model  using  known  methods. 
We  give  one  such  example. 

Proposition  2.1.  If  regular  interarrival  times  are  i.i.d.  (general  distribution  A)  with  rate  0  <  A  <  oo, 
and  service  times  are  i.i.d.  (general  distribution  G)  at  rate  p  >  0  and  are  independent  of  the  arrival  process 
then  { L(t )}  is  a  positive  recurrent  regenerative  process  if  and  only  if  X  <  p  where  p~l  d=  F{min(S,  Z)}  = 
Jo°°  -  G(s))ds. 

Proof:  From  Lemma  2.1  we  see  that  {■£.(*)}  is  stochastically  the  same  as  the  number  in  system  process  of  a 

FIFO  GI/GI/1  queue  with  inerarrival  time  distribution  A  and  service  time  distribution  that  of  min(Sn,Zn). 
The  result  thus  follows  from  well  known  results  on  such  queues  (see  for  example,  chapter  9  of  [10]).  ■ 

We  point  out  (more  generally)  that  if  ip  is  either  stationary  ergodic  or  governed  by  a  Harris  recurrent 
Markov  process  (HRMP)then  so  is  ip.  In  this  case  let  S  denotes  a  generic  t  .  time  from  a  Palm  version 
of  ip.  Then  {L(t)}  has  a  stationary  ergodic  version  as  long  as  A  <  £{min(S,  ~  2], [7], [8]). 

Our  second  model  assumes  that  the  superposition  of  negative  and  regular  arrivals  forms  a  renewal  process 
with  general  interarrival  time  distribution  F  assumed  to  have  non-zero  finite  first  moment;  T  will  denote 
a  generic  random  variable  ~  F.  Independent  of  the  past,  with  probability  0  <  p  <  1  an  arrival  is  regular 
and  with  probability  q  =  1  -  p  an  arrival  is  negative.  We  assume  that  customer  service  times  {5n  :  n  >  0} 
are  i.i.d.  ~  G,  where  G  is  a  general  distribution  with  non-zero  finite  first  moment,  p  d=f  (E(S))~1  is  the 
service  rate  (5  denotes  a  generic  service  time).  We  assume  that  the  service  time  sequence  is  independent 
of  the  exogenous  renewal  pivce®3  We  refer  to  this  as  th<*  rain  flipping  model  (CF).  Let  /V(t)  denote  the 
counting  process  of  a  non-delayed  renewal  process  having  cycle  lengths  i.i.d.  ~  G.  Let  71  denote  a  generic 
negative  interarrival  time.  T_  has  the  distribution  of  a  geometric  sum  (parameter  ?)  of  i.i.d.  F  distributed 
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interarrival  times;  E(T _)  =  E(T)/p.  E(N(T-))  =  /0°°  E(JV(t)){^)^Li  pn  1qF,n(di)}  denotes  the  expected 
number  of  potential  service  completions  during  71.  Here  7'"  denotes  the  n,h  -fold  convolution  of  F. 

Proposition  2.2.  For  the  coin  Sipping  RCS  model,  {£(*)}  is  positive  recurrent  regenerative  if  and  only  if 
6  =f  L-  E(N(T.))  -  1  <0. 

Proof:  First  we  prove  sufficiency.  Let  tn  denote  the  time  of  the  nth  negative  arrival.  Let  Xn  =  L(tn+).  It 

is  easily  seen  that  X  is  an  aperiodic  irreducible  Markov  chain  with  state  space  the  non-negative  integers.  We 
will  show  that  X  is  ergodic  by  use  of  Foster's  negative  drift  criterion  ( see  for  example,  Theorem  6.1  of  [9]).  To 
this  end,  it  suffices  to  show  that  for  some  k  there  exists  an  e  >  0  such  that  for  l  >  k,  Ei{X i)  <  l  —  c.  Observe 
that  £  is  precisely  the  expected  number  of  regular  arrivals  between  two  consecutive  negative  arrivals.  For  / 
sufficiently  large  we  thus  have  Et(X i)  -  I  as  S  <  0  by  our  hypothesis.  It  follows  that  for  any  0  <  c  <  |£  - 
E(N(T-))~  1|  a  k  can  be  found  giving  us  the  desired  negative  drift.  Thus  X  is  ergodic.  Assume  that  Xo  =  0, 
let  rn  denote  the  consecutive  times  at  which  X  hits  0  and  embedd  them  into  continuous  time  via  sn  =  tT%. 
By  ergodicity  E(n)  <  oo  and  hence  by  Wald’s  equation  E($i)  =  E(ti)E(T-)  =  E(rx)E{T)/p  <  oo.  Thus 
L(t)  is  positive  recurrent  regenerative  with  respect  to  the  renewal  process  (s„).  For  necessity,  we  will  apply 
(iii)  of  Theorem  11.3  of  [9].  Assume  that  8,  _  E(N(T-))  —  1  =  6  >  0.  We  mu>t  find  a  bounded  non¬ 
negative  function  g  such  that  Eig(X\)  >  g(l),  x  >  k  and  g(l)  >  sup{p(y)  :  y  <  k},  l  >  k.  For  k  large, 
we  will  use  g(x)  =  min(x,  Jfc  +  6).  Let  A  be  distributed  as  the  number  of  regular  arrivals  during  T_.  If 
/  >  k,  g(l )  =  min(/,  k  +  6)  and  Eig(X i)  =  E{min((l  +  A  -  N(T _)  —  1)+,  k  +  5)}  m  k  +  6  as  ib  — ♦  oo. 
Thus  Xn  is  transient.  In  particular  only  finitely  many  regular  arrivals  ever  find  an  empty  system.  But 
visits  of  X„  to  any  fixed  state  »  also  can  occur  only  finitely  many  times,  thus,  L{t)  — ►  oo  as  (  — •  oo 
and  L(t)  can  not  be  positive  recurrent  regenerative  with  respect  to  any  embedded  renewal  process.  For 
the  case  6  =  0,  we  can  apply  Theorem  8.1  of  [9]  with  y(i)  =  x,  A  —  [0,k]  and  d  =  p/q  to  show  that 
the  set  [0,*]  is  null  recurrent  for  all  sufficiently  large  k.  This  implies  that  X  is  null  recurrent  and  hence 
P(«i  <  oo)  —  1  and  E(sj)  =  oo  (via  Wald’s  equation).  In  this  case,  Let  p*  denote  the  limit  as  t  — ♦  oo  of 
1  jt f*P(L(s)  =  k)ds.  From  renewal  theory  p*  =  0  if  Ef^'  I(L(s)  =  k)ds  <  oo.  To  establish  this,  consider 
an  alternative  system  (designated  by  L  etc.)  with  the  same  input  except  each  service  time  has  been  scaled 
by  a  >  0.  For  cr  <  1,  L(t)  <  L(t)  for  all  t  Moreover,  for  a  sufficiently  small  6  <  0  so  that  L(t)  is  positive 
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recurrent.  Hence  EJq'  I(L(s)  <  k)ds  <  E /0*‘  I(L(s)  =  k)ds  <  oo.  Thus  p*  =  0  for  all  k  implying  that  L(t) 
is  not  stochastically  tight  and  hence  can  not  be  positive  recurrent  regenerative  with  respect  to  any  embedded 
renewal  process.  ■ 

Remark(2.1).  In  the  case  of  Poisson  regular  arrivals  (rate  A),  negative  interarrival  times 
i.i.d.  ~  A  (general  distribution),  and  i.i.d.  service  times  (general  distribution  G)  one  can 
use  the  method  used  in  proving  Proposition  2.2.  to  obtain  the  following  result:  L(t)  is 
positive  recurrent  regenerative  if  and  only  if  A E(T-)  —  E(N(T-))  —  1  <  0  where  now 
T_  ~  A. 

3.  Removal  of  the  customer  at  the  tail  of  the  queue 

In  this  section  we  consider  a  model  where  at  a  negative  arrival  epoch,  the  customer  m  the  system  who 
arrived  mo6t  recently  is  removed;  if  the  system  is  empty  then  nothing  is  done.  Observe  that  the  only  time  a 
customer  is  removed  from  service  is  when  he  is  the  only  one  in  the  system  at  a  negative  arrival  epoch.  This 
model  amounts  to  removing  the  customer  at  the  tail  of  the  queue;  hence,  we  refer  to  it  as  the  RCT  model. 
We  assume  that  the  superposition  of  regular  and  negative  arrival  times,  {f„},  forms  a  renewal  process  at 
rate  0  <  A  <  oo  with  interarrival  times  T„  i.i.d.  ~  F.  Let  Kn  =  1  if  the  n,h  arrival  is  positive  and  0 
if  negative.  We  shall  assume  that  {K„}  is  an  ergodic  Markov  chain  independent  of  {/„}  and  having  the 
property  that  P{K\  —  l|/f0  =  0)  =  1,  that  is,  two  or  more  negative  arrivals  cannot  occur  in  a  row.  Let 
p  =  P{Kx  —  l|/f0  =  1),  q  =  P(K\  =  0|Ao  =  1)  =  1  —  p;  0  <  p  <  1.  Service  times  S„  are  assumed  i.i.d.  ~  G, 
independent  of  {t„, /fn}  and  we  assume  that  0  <  p~l  d—  E(S )  <  oo. 

Proposition  3.1,  For  the  RCT  model,  { L(t )}  is  positive  recurrent  regenerative  if  and  only  if  A  < 

Proof:  For  mathematical  convenience,  we  imagine  that  every  arrival  (regular  or  negat  ive)  brings  a  service 

time  with  them;  Sn  denotes  the  service  time  of  the  nth  arrival.  Let  V(t)  denote  the  total  amount  of 
work  in  system  at  time  t.  Define  Wn  —  V(tn- )  if  K„  =  1;  Wn  =  V^(tn+)  if  Kn  =  0.  Define  An  = 
SnKn  -  T„  -  S„(l  -  K„+ 1).  It  follows  (with  a  little  thought)  that 

Wn+l  =  {Wn  +  An}+,  (3.1) 
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and  hence  we  have  a  reflected  random  walk  with  increments  A„  governed  by  the  Harris  ergodic  Markov 
chain  Xn  =  (Sn,T„,  Kn,  Kn+ 1).  From  known  results  (see  [7])  we  obtain  Harris  ergodicity  of  the  Markov 
chain  Zn  =  (Wn,Xn)  if  and  only  if  Ev{ A„)  <  0  where  v  is  the  stationary  distribution  for  X„.  It  is  easily 
shown  that  P„(Kn  =  0)  =  p/(l+p)  and  hence  (since  { Kn }  is  assumed  independent  of  {(Sn,rn)})  that 
F„(An)  =  — £(7*).  Thus  (W„,X„)  is  positive  Harris  recurrent  if  and  only  if  A  <  ;  in  this  case 

let  ic  be  the  stationary  distribution  for  (Wn,Xn).  P*{Wn  =  0)  >  0  and  hence  either  Pr{Wn  =  0,  K„  =  1)  >  0 
or  P*(Wn  =  Q,Kn  =  0)  >  0.  But  if  the  event  {W„  =  0 ,Kn  =  0}  occurs  (a  negative  arrival  leaves 
behind  an  empty  system)  then  so  will  { Wn+1  =  0,  Kn+\  =  1}  (the  next  arrival,  necessarily  regular,  finds 
the  system  empty)  since  by  assumption,  two  negative  arrivals  can  not  occur  in  a  row.  It  follows  that 
P*{Wn  =  0,K„  =  1)  >  0  and  hence  by  Wald’s  equation  L(t)  is  positive  recurrent  with  regeneration  points 
the  consecutive  times  at  which  a  regular  arrival  finds  an  empty  system.  Analogous  to  the  proof  of  Proposition 
2.2,  {L(t)J  is  not  stochastically  tight  when  E„{ A„)  >  0,  and  hence  can  not  be  positive  recurrent  regenerative 
with  respect  to  any  embedded  renewal  process.  ■ 

Remaik(3.1):  Results  like  Proposition  3.1  can  be  obtained  in  a  more  general  set-up:  when 
(Sn,T„,  Kn)  is  assumed  governed  by  a  positive  Harris  recurrent  Markov  chain  0  (and 
P(Kn+ 1  =  0| Kn  —  0)  =  0).  In  this  case,  however,  the  regenerative  structure  of  L(t)  is 
more  complicated. 

Remark(3.2):  When  the  input  to  RCT  is  a  general  marked  point  process  ip  =  (tn,  Sn,Kn), 
the  stability  conditions  can  be  quite  varied  as  the  following  example  shows.  Start  with  an 
initially  empty  M/M/ 1  queue  with  arrival  rate  A  and  service  rate  p;  p  =  A/p  <  1.  Let 
ip  =  (t„,Sn)  denote  the  corresponding  point  process  of  arrival  and  service  times.  Let  /„ 
denote  the  length  of  the  nth  idle  period  and  choose  U„  such  that  ~  U niform( 0,  /„) 

and  (L'n.fn)  are  i.i.d.  As  soon  as  the  n‘h  idle  period  begins,  let  a  negative  arrival  occur 
Un  time  units  later.  Construct  a  new  marked  point  process  ip  —  {Tn,^n,  Kn)  where 
denotes  the  nth  arrival  (regular  or  negative),  the  nth  arrival’s  service  time  (we  imagine 
that  negative  arrivals  bring  one  with  them),  and  K„  is  defined  as  before.  By  construction, 
every  negative  arrival  finds  an  empty  system  and  hence  has  no  effect  on  the  system;  p  <  1 
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remains  the  correct  stability  condition.  Moreover,  the  proportion  of  arrivals  that  are 
negative  is  strictly  greater  than  zero.  Finally  observe  that  the  point  process  regenerates 
at  each  negative  arrival  epoch. 
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